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ABSTRACT 


In this work, the performance of the NHPP-based software reliability model applying the Inverse-type 
distribution, which is widely utilized to various types of reliability life distributions, was newly identified. 
Accordingly, software failure time data was used to analyze reliability performance by predicting failures 
that may occur in the software operation, and the solution of parameters were estimated using maximum 
likelihood estimation. As a result, first, as a result of evaluating the criteria value (MSE and R7) for efficient 
model selection, the efficiency of the Inverse-Exponential model was evaluated as the best. Second, as a 
result of analyzing the attributes data (m(t), A(t), R(t)) that determine reliability performance, the Inverse- 
Exponential model was the most efficient. In conclusion, through various comparative analyses, the Inverse- 
Exponential model was found to have the best performance. Through the results of this study, the reliability 
performance of the Inverse-type life distribution for which there is no existing research data was newly 
analyzed, and basic design and test data necessary for an efficient software development process could also 
be presented. In the future, after exploring applicable statistical distributions for each software convergence 
industry, follow-up studies to find an optimal model will be needed. 
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Performance 


1. INTRODUCTION reliability analysis method to determine the 


goodness of fit after parameter estimation, Malik and 


In the current era of digital convergence, 
software technology is widely spreading and being 
utilized to various related industrial fields. 
Accordingly, the need for reliable software is rapidly 
increasing. Therefore, developers are currently 
concentrating on reliability studies to improve 
software quality. That is, the problem of improving 
the reliability of software becomes the most 
important issue for software developers. For this 
purpose, reliability models applying Non- 
homogeneous Poisson Process (NHPP) have been 
studied in various forms. Among these studies, the 
NHPP-based model applying the reliability 
performance attribute is attracting attention [1]. 
Also, regarding the Inverse-distribution proposed in 
this study, Pavlov and Lliev, Rahnev, Kyurkchiev 
[2] presented a method for calculating the error of 
the optimal approximation based on a modified 
Inverse-Exponential software reliability model. 
Fatima and Ahmad [3] proposed an improved 
Inverse-Rayleigh distribution by applying a new 


Ahmad [4] proposed an improved model of the 
Inverse-Rayleigh distribution using Alpha Power 
Transformation. Therefore, Voda [5] explained with 
an example that the Inverse-Rayleigh distribution is 
applicable to various lifetime distributions. Also, 
with respect to the software reliability model, Prasad 
and Rao [6] analyzed the performance of the NHPP 
model applying the Inverse-Rayleigh distribution, 
and Huang [7] evaluated the NHPP software 
reliability properties by applying the performance 
attribute function. Kim [8] compared the efficiency 
in terms of statistical process control for the 
reliability attributes of the NHPP model using 
Inverse-Rayleigh and Rayleigh distributions, Yang 
[9] also defined the performance of NHPP software 
reliability models using Exponential-type 
distributions. Additionally, many researchers are 
working on exploring the best distribution by 
applying various types of lifetime distributions to 
studies related to software reliability. 
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Accordingly, this study presented a new 
NHPP-based software reliability model using the 
Inverse-type life distribution, which has been proven 
to be efficient in reliability testing of various life 
distributions. Along with this, the reliability 
attributes of the proposed model were explored and 
its performance was newly identified. Also, we will 
propose an optimal model through the analyzed data. 


2. RELATED RESEARCH 
2.1.1 NHPP model 


The NHPP is known as a probability-based 
distribution model that predicts future failures based 
on the number of failures that occur at a given area. 
Thus, this study aims to evaluate the performance of 
the reliability model using this probability-based 
NHPP model. That is, the NHPP model has been 
widely used to model the number of failures N(t) 
found between observation times (0, t) in reliability 
measurement. In a software system where failure 
times occur at different intervals and failures occur 
continuously, if the number of failures occurring per 
unit time is N(t), then N(t) follows the Poisson 
distribution and also satisfies inhomogeneity. 


Accordingly, the NHPP model can be said to 
be a probability-based model that can predict 
software reliability based on the number of failure 
occurrence. Thus, using these properties of the 
NHPP model, it can be defined as follows. 


[mn]? 2e-" © 
n! 
Note that n = 0,1,2,-:+ ©. 


P{N(t) =n} = (1) 


where m(t) refers to a mean value function that has 
the property of estimating the true value and can be 
defined as Equation (2). Thus, the intensity function 


X(t), which has properties representing the 
instantaneous failure occurrence rate, can be 
developed as follows. 
t 
m(t) = A(s)ds (2) 
0 
dm(t) 
=A(t 3 
ie 740 (3) 


2.1.2 NHPP software reliability model 


In this work, we seek to solve the cost 
attribute problem of the proposed NHPP software 
development model based on failure time collected 
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during normal operation. This study reflects the 
failure phenomenon of generally developed software 
and aims to study it based on finite failure, in which 
no further failures occur after the failure is repaired. 


Therefore, this study is intended to be 
developed based on the finite failure NHPP model 
by reflecting realistic failure situations. Accordingly, 
applying Equations (2) and (3), the attribute 
functions representing the performance of the cost 
model are as follows [10]. 


m(t|@,b) = OF (t) (4) 
A(t]O,b) = OF (t)’ = Of (t) (5) 


Note that 0 is the residual failure rate, and F(t ) is the 
cumulative distribution function. 


From the result derived above, m(t) represents the 
performance that can predict the true value, and A(t) 
represents an attribute that represents the intensity at 
which a failure may occur. 


Accordingly, the likelihood function of the NHPP 
model applying the attribute functions m(t) and 
A(t) can be developed as follows. 


Lyupp(0|x) = (T] ro) exp|—m(xy,)] (6) 


Note that x = (1, X2,X3°** Xn) 


2.2 NHPP Goel-Okumoto Basic Model 


The Goel-Okumoto model is widely known 
as the most basic NHPP model because it is based on 
the basic concept that the number of faults 
discovered per unit time is proportional to the 
number of faults remaining at that time. Also, 
because the time distribution for failure occurrence 
per defect in the Goel-Okumoto basic model has 
exponential distribution characteristics, it is also 
called an exponential-type basic distribution model. 


Therefore, if the expected value of the defect causing 
the failure in the finite failure situation of this model 
is expressed as 0 and the defect search rate is b, it can 
be developed by considering b as a fixed constant. 


Since the failure rate can be considered a constant 
with a certain form, the performance function 1s as 


follows [11]. 


m(t|6, b) = OF(t) = @(1 —e-"*) (7) 
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A(t|0,b) = Of (t) = Obe~” (8) 


Accordingly, the log-likelihood function of this 
distribution model calculated by applying Equation 
(6) 1s as follows. 


n 
InLyupp(O|x) = nlnd +nInb —b ». +i 
k=1 
—@(1 —e-P*n) (9) 


Finally, the parameter estimator (@y,¢, byzz) of the 
Exponential-basic model to be obtained in this work 
can be calculated using maximum likelihood 
estimation (MLE) to Equation (9) and then using the 
bisection method. Therefore, Equations (10) and 
(11) show the final calculation equations for 
calculating the parameters. 


OlnLyypp(O|x) _ 


= 2-14.75 = 9 (10) 
00 6 


dlnL O|x nw 2 . 
cere (OND =e > Xn — Ox,e-?* =0 (11) 


i=1 


2.3 NHPP Inverse-Exponential Model 


The Inverse-exponential distribution has a 
bathtub-shaped risk rate function, so it is a useful 
distribution for the analysis of reliability life data 
with the characteristic that the risk rate varies with 
time. In particular, the Inverse-exponential 
distribution is known as a distribution suitable for 
load-intensity reliability, which represents the 
probability that the system will operate normally 
when stress is applied to the system stochastically. 
Accordingly, the Inverse-exponential distribution 
plays a very important role in measuring the 
reliability of a system in the field of reliability. 


Therefore, the function F(t) can be defined as 
follows. 


F(t) = e~@0* (12) 


If the functions obtained above are substituted into 
Equations (4) and (5), the performance attribute 


functions of this model are as follows [12]. 


m(t|6,b) = @e~@2)* 
A(t|O, b) = Ob-1t-Ze-&* 


Thus, the log-likelihood function of this NHPP 
model calculated by applying Equation (6) is as 
follows. 


InLyypp(O|x) = nlné — nlnb (15) 


n n 
+2 > xj - >. ox" — Ge-(bxn)* = 0 
i=1 i=1 


Therefore, the parameter estimator (OQy,2, Dupe) of 
the Inverse-exponential model can be calculated by 
applying MLE to Equation (15) and then using the 
bisection method. Accordingly, Equations (16) and 
(17) show the final calculation equations for 
calculating the parameters. 


OlnLyypp(O|x) n eon ae 
Se —(bxn) = 0 16 
00 a ° — 
n 
alnL O|x n 1 
ener OO = 2 — (17) 
b b = 
1 | 
—6 -(x%n)~" = 9 
b2x,, 


2.4 NHPP Inverse-Rayleigh Model 


Like the Rayleigh distribution, which is 
known to be a suitable model in the field of system 
lifetime testing, the Inverse-Rayleigh distribution is 
also a life distribution that has many applications in 
the field of software reliability. In particular, the 
Inverse-Rayleigh distribution is a model that has 
been proven to be efficient in reliability analysis of 
various life distributions and has been confirmed to 
be suitable for software reliability testing. 


Accordingly, after applying these characteristics to 
reliability research, many researchers confirmed that 
this model can be used as a life distribution in 
reliability test and property analysis as follows. 


F(t) = exp (- =) (18) 


f(t) = S exp (- =) (19) 


Therefore, if the functions obtained above are 
substituted into Equations (4) and (5), the 
performance attribute functions of this model are as 
follows [13]. 


Journal of Theoretical and Applied Information Technology 
15" May 2024. Vol.102. No 9 


SZ 


© Little Lion Scientific 


ISSN: 1992-8645 


b 
m(t|6,b) = F(t) = dexp (- =) (20) 


2 


b b 
A(t|O,b) = f(t) = @ | = exp (- -.) | (21) 
Therefore, if arranged in the same way as Equation 
(15), the log-likelihood function of this model can 
be written as follows. 


InLyypp(O|x) = nln2 + nlné + nlnb 


“4 4 b 
+b) In 3 -b ) =~ bexp(-5) (22) 
i=1 . i=1 ue 


That is, the parameter estimator (Oy,2, Dyzz) of the 
Rayleigh model to be obtained in this work can be 
calculated by applying MLE to Equation (22) and 
then using the bisection method. Thus, Equations 
(23) and (24) show the final calculation equations for 
calculating the parameters. 


AlnLyypp(O|x) n b 
70 eas 2 (23) 
n n 
OlnLyypp(O|x) n 1 1 
r= ==+) In(s}-) 5 
db b ¢ xp) Lax} 
t=1 t=1 
6 b 
+ —, exp “ee = 0 (24) 


3. RELIABILITY PERFORMANCE 
ANALYSIS 


In this work, the performance properties 
applying Exponential-type life distribution model 
were analyzed by the step-by-step sequence of the 
presented solution as follows. Also, the optimal 
model was presented based on the analyzing data. 


Table | [14] shows the software failure time 
data cited in this paper. This data refers to the 
collection of failure times that occurred while 
operating the software system. 


Also, this data is a collection of 30 failures for a total 
of 187.35 hours, which occurred due to design and 
analysis errors in the software development process. 
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In this study, Laplace trend analysis was used 
to judge whether the software failure time cited were 
applicable to this work. 


Table 1: Software Failure Time Data. 


Failure time Failure time 
(hours) (hours)x 107? 


Failure 
number 


= 
=~ Soe marniant Wn = 


NNN NNN NN WN N YF RR FPR FF Re 
OMANI NHN PWN RK TO WON HW Nn FW WN 


eS) 
a) 


In general, if the Laplace trend analysis result of the 
cited data is distributed between '-2 and 2', this data 
is said to be reliable. 
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Figure | shows the results of the Laplace trend test 
analyzed by applying the data presented in Table 1. 
That is, it can be seen that all result data exists 
between '0 and 2'. Accordingly, it can be said that the 
cited software downtime data is applicable to this 
study. 


Laplace trend test 


agi: a 

$ 9 ve 
0.24: ; . ® 
@ 2 


59; 9 Ib 13 15 17 19 21 2 2 27 2 


Laplace factor 


0” 
a 
Failure time 


Figure 1: Results of Laplace Trend Test. 


The parameters (8, b) of the NHPP model were 
calculated using the MLE as shown in Table 2 [15]. 
Therefore, among the parameters of the applied 
model, 6 is a software residual failure, and b is a 
shape parameter that makes the shape of the 
applied life distribution. 


In this study, we will also analyze the R* and MSE, 
which are criteria for determining an efficient model. 


Table 2: Parameter Solution Using MLE. 


NHPP 
model 


Type 


Goel- 
Oia 32.9261 | 0.1297 


Inverse- 
Inverse-type Exponential 41.2881] 0.1692 
life 
Inverse- 
Rayleigh 30.0100] 1.6520 


R*, which is used as a standard for explaining the 
difference between actual values and observed 
values, is expressed as follows. 


Basic model 


distribution 
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> _ (mee) - me)? 


R= 1 Sy 
D,,,., mi) ~ Djs ™)/n)) 


(25) 


Note that M(x;) is the cumulative number of failures 
estimated from m(t). 


In comparison, if the coefficient of determination is 
large, the error is small and it is considered a 
relatively useful model. In other words, eventually 
the error becomes smaller, so it is considered a 
relatively efficient model. 


MSE is a standard for comparing the difference 
between the real value (actually observed value) and 
estimated value (predicted value) and is as follows. 


>, _ (med - a)! 


MSE ==) —__ (26) 


Note that n used in this equation is the number of 
observed failures. 


Figure 2 is the result of analyzing the model 
properties using MSE, and this study intends to use 
this value as reference data to determine the 
suitability together with the efficiency. 


Mean Square Error Vs. Failure number 


«eee Goel-Okumoto quae Inverse-Exponential 


Inverse-Rayleigh 


A ‘ ‘nt See oe ae ee ee 


Se og go, ga 
SF xt Ft FH HF NNN NN ON 


Failure number(times) 


Figure 2: Analysis of MSE. 


When selecting an efficient model, the smaller the 
value of MSE, the smaller the error predicting the 
true value, so it is determined as a relatively 
efficient. 
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Table 3 shows detailed data analyzing the change in 
MSE value with failure time to verify the efficiency 
of the model [16]. 


In general, if the coefficient of determination 
is greater than 0.8 (80%), this model is said to be 
efficient. As shown in Table 4, among the models 


Table 3: Analysis Data Using MSE. 


Failure MSE 
number 
Okumoto | Exponential Rayleigh 
4 0.1555 0.3276 4.6671 
5 0.7089 0.0180 12.2861 
1.7288 14.8474 
2d 0.2066 0.1894 0.2887 
28 0.0877 0.0815 0.1222 


proposed in this work, the Inverse-exponential and 
Goel-Okumoto models are judged to be efficient. In 
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other words, the Inverse-Exponential model is 
judged to be more efficient and useful than the 
Inverse-Rayleigh model. 


The m(t), which refers to the reliability 
performance properties of the proposed model, is an 
important function to measure _ reliability 
performance. In particular, the m(t) function 


Table 4: Model Efficiency. 
NHPP 2 
Basic model | Coch 0.8956 | 32.9379 
Okumoto 
0.9359 | 20.2035 
Exponential 
Inverse- 
Ravin 0.4747 | 165.750 
represents the expected value of software failure 


occurrence and is also an important indicator of the 
predictive power of estimating the true value. 


Inverse-type 
lifetime 
distribution 


Table 5 is a simplified summary of the equations for 
calculating the m(t) [17]. 


Figure 3 shows the trend of prediction ability 
to estimate the true value over the passage of failure 
time. 


Table 5: Mean Value Function (m(t)). 


6(1 —e~"*) 


Type NHPP 
yP model 
Goel- 


Basic model Giaiats 


Inverse- 


—(bt)"* 
Exponential me 


Inverse-type 
lifetime 
distribution | Inverse- 


Rayleigh 


bexp (- 7) 


When analyzing the trend curve, all models show 
results that do not accurately predict the true value 
but estimate the error value. However, in analyzing 
the performance of predicting the true value, the 


3987 


Journal of Theoretical and Applied Information Technology 
15" May 2024. Vol.102. No 9 


SZ 


© Little Lion Scientific Ga 


ISSN: 1992-8645 


Inverse-Exponential model with the smallest error 
can be said to be the most efficient. 


Mean value function Vs. Failure time 


«eee Goel-Okumoto eum Inverse-Ex ponential 


oem Inverse-Rayleigh oes True Value 


30 
25 + 
20 + 
my 
— 45 
E 
10 - 
5 
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xn | WN * mon mm vw ~~ Ww 
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Failure time(hours) 


Figure 3: Performance Analysis of m(t). 


The A(t), along with the m(t), is an important 
function to measure the reliability performance 
properties. In particular, the A(t) is a failure rate 
function, which means a failure rate per fault and is 
also an important index indicating the intensity of 
software failures. 


Therefore, Table 6 briefly summarizes the equations 
for calculating the intensity function [18]. 


Figure 4 shows the trend of failure rate 


Table 6: Intensity Function (A(t)). 
NHPP 
Basic Goel- =e 


Inverse- 

Inverse-type | Exponential 
lifetime 
distribution | Inverse- 


Rayleigh 


occurrence intensity over the entire failure time 
range. Therefore, as a result of analyzing the failure 
rate trend of the proposed NHPP models, similar to 
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general failure phenomena, the intensity function 
initially increased, but as time elapsed, the failure 
rate was removed and the intensity function showed 
an efficient trend with a gradually decreasing 
pattern. 


That is, the Inverse-Exponential model showed the 
lowest failure rate and was very efficient, but Goel- 
Okumoto model showed inefficiency that only 
decreased continuously. 


Intensity function Vs. Failure time 


eee Goel-Okumoto eee [nverse-Exponential 


ome Inverse Rayleigh 


20 + 
16 + 
=. Qe. 
= 
oa 
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Failure time(hours) 


Figure 4: Performance Analysis of A(t). 


Table 7 shows data values analyzed in detail 
according to the number of failures that occurred 30 
times using attribute functions (m(t), A(t)) that 
represent reliability performance, which is the core 
topic of this study. 
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Table 7: Trend Analysis Values of Reliability Performance Attributes. 


Reliability Performance Attributes 
Failure 


x 
(hours) 


ULT9 25.78026724 24.99860162 | 29.65479731 0.926814509 1.064874711 0.059952836 
12.536 26.44852342 25.76762543 29.696 18234 0.840141683 0.969072062 0.049803941 


12.973 26.80545455 26.18011141 29.7168665 0.793847715 0.91936993 1 0.044969887 


28.34272074 27.98935819 | 29.79626959 0.59446429 0.715704899 0.028016511 


20,595276/3 28.29503943 29.80800676 0.561707779 0.683654139 0.025743212 


15.98 28.78210766 28.52345327 29.8 164833 0.537475806 0.660158462 0.024141605 
16.385 28.99416702 28.78540434 | 29.82590261 0.509971707 0.633693254 0.022402387 


29.276/73326 29.13958574 | 29.83813881 0.473322866 0.598730413 0.020208507 


29.40551645 29.30322664 | 29.84360286 0.456619686 0.582896868 0.019253347 


29.567428 29.51118859 | 29.85037786 0.435619758 0.563068208 0.018090554 
18.122 29.78729468 29.79802889 29.8594184 0.40710305 0.53625941 0.016576888 
18.735 30.02718608 30.11770269 | 29.86908831 0.375989135 0.507123911 0.015007243 
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The reliability function ( R(t)) is an 
important function that can predict future reliability 
performance along with the attribute function (m(t), 
Mt)) analyzed above. In particular, the purpose of the 
reliability function is to analyze the trend of future 
reliability by assigning a random duty time again 


E-ISSN: 1817-3195 


Table 8 shows the result of analyzing the 
reliability performance trend in detail after putting 
future mission time into the NHPP models proposed 
in this study. For reference, the mission time 
( 145Hx107' ) presented in Table.5 was 
numerically converted to facilitate calculation. 


after the final failure time (x,, = 18.735). 
Therefore, in this study, after putting mission time TapICS Ana eae: iG): 
into the proposed model, we try to predict and 
evaluate the future reliability performance. Here, 
reliability means the probability that a failure occurs 
at the test point and no failure occurs between the 
confidence intervals. Therefore, future reliability 
(R(t)) can be defined as follows [19]. 


Mission Reliability Function R(t) 
Time 


Exponential | Rayleigh 


0.496273344 | 0.985038791 | 0.974507874 
0.43125056 | 0.977317633 | 0.969311894 
0.378066466 | 0.970006765 | 0,964494471 


RG) = exp[—{m(x, a T) — m(X,)5] 
= exp[—{m(18.735 + tT) —m(18.735)}] (27) 


Note that t is the mission time. 


As shown in Figure 5, as a result of analyzing the 
reliability trend after putting in the mission time, the 
Goel-Okumoto model can be said to be inefficient 
because the reliability decreases as time goes by. 
But, the Inverse-Exponential and Inverse-Rayleigh 
model, which show a consistently high and stable 
trend compared to the Goel-Okumoto model, can be 
defined to be very efficient. 


0.184611493 | 0.928013092 | 0.938792645 


0.17070028 0.923067953 0.935991589 
0.15861571 0.91833557 0.93335656 


0.14806775 0.913802561 0.930874571 
eS) 0.138819776 | 0.909456638 0.928533889 
9 0.13067754 0.905286501 0.926323895 


Reliability Vs. Mission time 


«eee Goel-Okumoto eum Inverse- Exponential 


eww Inverse-Rayleigh 


oa 


reat Maggs 
01 4 Blan Oe 
0 serene eee 

con noene ne eros s 


a - = oF 
Mission dicaiatbameaal 

13.0 0.087693741 | 0.874181279 0.910976334 
14 0.084965794 | 0.871307326 0.909661737 


14.5 0.082486226 | 0.868528654 | 0.908407738 


Ss IDs SE OSDn Oh SR A A SW | 


Figure 5: Performance Analysis of R(t). 
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The topic of this study was the analysis results of 
model efficiency (MSE, R*) and attribute data (m(t), 
X(t), R(t)) that have a significant impact on model 
performance. Also, Table 9 shows the final 
evaluation results after comprehensively comparing 
the performance attribute data of the proposed model 
based on the research data developed in this work. 


Accordingly, if this research data can be 
utilized efficiently in the early stages of software 
development, it is believed that this data can be 
helpful not only as basic design data needed by 
developers but also as attribute data required to 
improve reliability [20]. 


Table 9: Reliability Performance Evaluation. 


Bi Performance Attributes 
NHPP iciency 
model R 


Goel- 
Okumoto 


Good | Worst 


Inverse- 
Exponential 


TIWETSE Worst | Worst | Worst | Worst | Best 
Rayleigh 


4. CONCLUSION 


If a software developer can design a 
reliability prediction model with failure time data 
collected in the early stage of analyzing and testing 
a program, developers will be able to predict 
software failure time in advance to increase 
reliability and ultimately improve software quality. 
Thus, the performance of the NHPP reliability model 
applying Inverse-type life distribution property, 
which has been widely known to be suitable for 
reliability analysis, was analyzed and its attributes 
were identified. 


The results of this work are as follows. 

First, as a result of analyzing the reference data 
(MSE and R*) for efficient model selection, the 
efficiency of the Inverse-Exponential model was 
evaluated as the best. 


Second, as a result of analyzing the performance 
attribute data (m(t), A(t)), the Inverse-Exponential 
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model with excellent predictive ability of true value 
and low failure rate was the most efficient. 


Third, as a result of the reliability test, the Inverse- 
Exponential and the Inverse-Rayleigh model, which 
showed consistently high and stable reliability, were 
efficient. However, the Goel-Okumoto model, 
which showed the attribute of continuously 
decreasing reliability with mission time, was 
inefficient. 


In conclusion, this study can present solution 
techniques and basic design data that can analyze 
and predict performance attribute data needed by 
developers during the early software development 
process. Additionally, follow-up research will be 
needed to use the results of this study to find an 
optimized reliability model suitable for related 
software industry fields and to explore attribute data 
related to reliability performance. 
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